A1 = readmatrix('NiI2_MH.dat')

x1 = A1(1:281,4)
y1 = A1(1:281,59)*130


hold on
%scatter(x1,y1,5,[0 0 0.8],'filled')
plot(x1,y1,Color=[0 0 0])
xlabel('Magnetic Field [Oe]')
ylabel('chi [emu/mol]')
ylim([-0.06,0.06])
xlim([-80000,80000])
yticks([-0.04:0.02:0.04])
xticks([-60000:20000:60000])

set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 300, 210])